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1. Introduction 

Spontaneous particle production in the presence of strong electric fields is perhaps one 
of the most interesting phenomena associated with the charged quantum vacuum [H[2]. 
The study of supercritical effects and resonant particle production by strong Coulomb- 
like potentials dates back to the pioneering works of Pieper and Greiner [3] and Gershtein 
and Zeldovich [I] where it was shown that spontaneous positron production was possible 
when two heavy bare nuclei with total charge larger than some critical value Z c collided 
with each other. The critical Z c is the value for which the 15* state of the hydrogenlike 
atom with total charge Z has energy E = —m. The idea behind supercriticality lies 
in spontaneous positron emission induced by the presence of very strong attractive 
vector potentials. The presence of a strong electric field induces the energy level of an 
unoccupied bound state to sink into the negative energy continuum, i.e, an electron of 
the Dirac sea is trapped by the potential, leaving a positron that escapes to infinity. The 
electric field responsible for supercritical effects should be stronger than 2m e c 2 , which 
is the value of the gap between the negative and positive energy continua. Such strong 
electric fields could be produced in heavy-ion collisions [U EJ [6] 

In the last years the one-dimensional Dirac equation in an external potential V has 
been studied in connection to the Levinson theorem in relativistic quantum mechanics 
[3, El E] - The interest in this problem can be found in their applications in semiconductor 
physics, field theory and one-dimensional QED among others. One-dimensional Dirac 
electrons can be considered as three-dimensional problems when V depends only on one 
space variable [7J. The relation between supercritical states and the Levinson theorem 
in the one-dimensional Dirac has been established by Ma et al [5], and Calogeracos and 
Dombey [8]. 

Recently the relation between transmission resonances and supercriticality effects 
have been studied in a one-dimensional barrier [10J and in a Woods-Saxon potential 
[TT| [12] , which is a smoothed form of a square well. Here the Dirac equation presents half- 
bound states with the same asymptotic behavior of those obtained with the potential 
barrier. Supercritical states with E = —m have been also studied in a one- dimensional 
cusp potential [131 E] an d in a class of short-range potentials [15] . 

In the present article we discuss the problem of computing the phenomenon of 
resonant particle creation in a one-dimensional symmetric cusp potential of the form 



where V > 0, a > and e is the charge of the particle. 

The potential ([1]) is an asymptotically vanishing potential for large values of the 
space variable x. The parameter V shows the depth of the well. The constant a 
determines the shape of the potential. The expression ([1]) can be regarded as a screened 
one dimensional Coulomb potential [16]. 

It is well known that a vectorial delta potential is strong enough to pull the bound 
state into the negative energy continuum E = —m [I7J[T8], nevertheless this supercritical 
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state does not evolve to a real resonant state. It is of interest to investigate the resonant 
behavior of the lowest energy states in a potential exhibiting a delta point interaction 
as asymptotic limit. The cusp potential (CD) reduces to an attractive vectorial delta 
interaction of strength g v in the limit 2Va — > Qui cLS CL — > 0. The attractive potential (pQ) 
does not possess compact support and, in contrast to the delta potential case, it can 
exhibit more than one bound state. 

It is the purpose of the present article to show that the cusp potential Q is strong 
enough to sink the lowest electron bound state into the negative continuum to create a 
resonant state. In order to calculate the process of embedding of the bound state into 
the negative energy continuum we construct the Jost solutions and their corresponding 
Jost functions [191 120] • We estimate, with the help of the phase shift and the Wigner 
time-delay [22l l2Tj . the mean life of the resonance and how it depends on the shape and 
strength of the potential. 

The article is structured as follows. In section 2 we solve the 1+1 Dirac equation in 
the potential (OQ) and derive the equation governing bound states, energy resonances and 
supercritical states. In section 3 we compute, with the help of the Jost functions, the 
transmission and reflection amplitudes and the phase shifts. We derive the condition 
for energy resonances and estimate the mean life of the antiparticle state. Finally, in 
section 4 we summarize our results. 



2. Solution of the Dirac equation and energy resonances 

In two dimensions, the Dirac equation in the presence of the spatially dependent electric 
field ([[]) can be written as 



{7°(<9 t - iV exp(- \x\ /a)) + ^d x + m] ^ = (2) 

where the 7^ matrices satisfy the commutation relation {7^,7^} = 2?7 / " / and the metric 
has the signature (—,+). Here and along the paper, we adopt the natural units 
h = c=l. 

Since we are working in 1+1 dimensions, it is possible to choose the following 
representation of the Dirac matrices 

Taking into account that the potential (CQ) does not depend on time, we can write the 
spinor \1> as follows 

Substituting the gamma 7^ matrices ([3]) into the Dirac equation (T2]), we obtain the 
following system of coupled differential equations 

+ i(Vexp(- \x\ /a) + E) J #1 + m^ 2 = (5) 
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The spinor solution of the system (jS!)-®, exhibiting a regular behavior as x — > +00, 
can be expressed in terms of the Whittaker function M k ^(y) [23] as follows 

where b\ is a constant, y = —2iaVexp(—x/a) and 



k = iEa- 1/2, y, = iaVE 2 - m 2 (8) 

For negative values of x, we have that the solution of the system ([S])-® showing a 
regular behavior as x — > —00 can be expressed in terms of the Whittaker functions 

M k ,n(y) as 

where C\ is a constant and y = — 2aiV exp(x / a) 

Imposing the continuity of the spinor solution ^ at x = 0, and demanding the 
existence of non-trivial values for c\ and 61, we obtain that resonances and bound levels 
satisfy the energy condition: 

Q + fi + k^j Ml +l41 (-2ieaV) + m 2 a 2 Ml^-lieaV) = (10) 



Given the values of the strength V and the shape parameter a, the expression (1T0I) 
exhibits solutions with E real, showing that the cusp potential ([1]) is able to bind 
particles, a property that, for instance, is absent in the pure one-dimensional Coulomb 
potential. The number of energy levels in the cusp potential decreases as the shape 
parameter decreases. 

As the potential depth increases for large values of V, the energy eigenvalues of 
the bound states decrease. When the lowest energy bound state reaches the value 
E = —m, this level merges with the negative energy continuum and the potential 
becomes supercritical. For E = —m, we have that y — and, using Eq. fflQj) . we obtain 
the condition for supercriticality 

M kfi (-2ieaV c )M ktQ (-2ieaV e ) - M k+li0 (-2ieaV c )M k+lt0 (-2ieaV c ) = 0(11) 

As soon as the potential strength V becomes larger than the supercritical value V c , 
the lowest bound state dives into the energy continuum becoming a resonance. Fig. [T] 
shows the dependence of the supercritical potential V c on the shape a. As a — > we 
have that — 2aV c approaches tt, showing that the Dirac equation in the presence of a 
vectorial delta interaction of strength g v becomes supercritical as g v — tt [17]. 

Fig. [2] shows, for a given value of the shape parameter a, the dependence of the 
resonant energy of the lowest bound state on the potential strength V. The curve of the 
imaginary part as a function of the real part of the resonant energy shows a concavity 
that decreases as the potential V increases. 
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Figure 1. Plot of the supercritical potential strength V (E = —m) against the shape 
a. 
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Figure 2. Resonant energy spectrum for a = 1.57 with 3.145533 < V < 3.15 
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3. Transmission amplitude and phase shifts 

In order to study resonant states associated with Dirac particles in a cusp potential, we 
proceed to discuss the scattering process with the help of the Jost functions [19, 22]. 
Resonances can be identified as poles of the scattering matrix S [T9(, [2TJ |22] • Since the 
components of the scattering matrix as well as the reflection and transmission coefficients 
R and T can be expressed in terms of the reflection and transmission amplitudes, we 
proceed to compute r and t with the help of the incoming ± (x) and outgoing ip Jost 
solutions, where plus and minus indicate the asymptotic right and left solutions. 

Since the cusp potential (0Q) decays faster than 1/x 2 for large values of x, there 
are solutions of the Dirac equation (jSJ) that behave asymptotically as free traveling 
plane wave spinors. Those solutions exhibiting this asymptotic behavior are called Jost 
solutions [22]. The one- dimensional Dirac equation possesses two outgoing ^ and two 
incoming ^ Jost solutions [201 121] 

Since the Jost solutions i/j±(E,x) and ip±(E,x), are linearly independent, we can 
write regular solutions ty(~> (E,x), (E,x) of the Dirac equation as [20] 



*W (E, x) = - —J— [W (E, x) , V± (E, x)) v± (E, x) = (12) 

-W (^ (±) (E, x) , V± (E, x)) ^± (E, x)] 

= - w(l 1 UF , f/| ±} (E) <p± (E, x) - (E) j, ± (E, x) 
W (V>±, <p±) [E) L 

where the Wronskian W (ift(x), <f(x)) of the two two-component solutions ip(x) and <p(x) 
of the Dirac equation ip(x) = (ipx(x),ip2(%)) and ip(x) = ((fi(x), ^(x))* is defined as 
[26JE5] 

W^(x)^(x))^det( ^ l( . X \ ^ ) = Mx)Mx) ~ Mx)Mx) (13) 

V W2\X) ( P2{X) J 

It is easy to see that if W(ip, (p) ^ we have that ip(x) and <p(x) are linearly independent, 
/_ (E), (E) are in Eq (TT2"j) , analogous to the non-relativistic case, the Jost 
functions [22]. 

We make use of the Jost solutions in order to describe a scattering process. The 
solution ip + should be a linear combination of the solutions ip- an d <£-, that is, 

t±(E) ip±(E, x) = <p T (E, x) + r T (E) ip T (E, x) (14) 

Since the cusp potential is even, using the invariance under parity [1] and the 
representation ([3]) for the gamma matrices we have that 

ip±(E, x) = -ia y ip T (E, -x), (f±(E, x) = -ia y ip T (E, -x) (15) 

The parity conditions (j!5p imply that the left and right transmission and reflection 
amplitudes are identical. The transmission amplitude t(E) = t^{E) can be expressed 
as 

M = *w = -Wf¥rm (16) 
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analogously, we have that the reflection amplitude r(E) = r T (E) reads 

r(E).r {E)- w{ ^ +){E y (17) 

The incoming Jost solution for the one dimensional Dirac equation ([2j) takes the 
form 



m 



(-2iaeVyj2E(E - V^ 2 - m 2 ) V f ^fc+i, M (2/J, 



a spinor that exhibits the following asymptotic behavior 

m ( \+ k -» \ 

(p~(x) — > ma exp(ixv E 2 — m 2 .) (19) 

^2E(E-VE 2 -m 2 ) \ 1 / 

We can see that the solution f|T8|) behaves asymptotically as a traveling to the right 
plane wave solution to the free Dirac equation. The outgoing Jost solution ip_ can be 
written as 

My) = = nM J M \ V 1 (20) 

(-2iaeV)-vJ2E(E + VE 2 - m 2 ) \ V IM h+i,-^[V), 



such a spinor exhibits the following asymptotic behavior as x — * — oo 



ip-(x) — > — - t — ( ™ a ] exp(— ix\/ E 2 — m 2 ). (21) 



2E{E + V^ 2 ^ 



This asymptotic form shows that the Jost solution (1201) asymptotically behaves like a 
traveling left solution wave to the free Dirac equation. The right outgoing Jost function 
i)+ is 



rn 



My) = , ™ % J 7, I ' ( 22 ) 

{-2iaeV)-»J2E{E- VE 2 -m 2 ) \ V ' %c IW 



such a spinor exhibits the following asymptotic behavior. 



if)+(x) — > — ; = ( ™ a ) exp(ixV-E 2 — m 2 ) (23) 



'2E(E - VE 2 -m r 
Finally, we have that the incoming Jost solution (p + takes the form 

{-2iaeVY\ 2E{E + V^ 2 - m 2 ) V ^ 



whose asymptotic behavior is that of an incoming wave from the right. 

m ( l+^+ fc 

tp+(x) — > — - t — ma ) exp(— ixVE 2 — m 2 ) (25) 

2E(E + VE 2 - m 2 ) V 1 
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Figure 3. Phase shift 6 as function of the real part of the energy E for V = 3.15, 
a = 1.57. Notice that the phase shift 6 of the resonance value %IE = —1.0032 is larger 
than 7r/2 



With the help of the Jost functions and using Eq. (fl6|) . we have that the 
transmission amplitude t for the Dirac equation (J2J) takes the form 

ma 2^M k+1 ^{-2iaeV)M k ^{-2iaeV) - ^^M k+l ^{-2iaeV)M k ^{-2iaeV) 

' > = 1 1 h i 71 h- \ 2 

5 + * (-2iaeV)^(Ml _„(-2iaeV) + il^Af^C-HaeV)) 

Analogously, we have that the reflection amplitude r(i£) is 

r (£) = j ma M fc+liM (-2mey)M fc+1 ,_ M (-2mey) - M fc| _ M (-2mey)M M (-2mey) 

* + k + * (-2zaeV)^(Ml^(-2 ia eV) + &^Ml +1 ^{-2iaeV)) 
It is not difficult to see that the components t(E) and r(E) of the scattering matrix 
S have the same poles, they are the roots of the Jost function f+(E) 

f + (E) =W(1>-,il> + ) (E) = Ml_^-2iaeV) + ^ + ^ 2 ~^ M 2 k+1) _^(-2iaeV) (28) 

It is worth noticing that the equation f+(E) = coincides with the resonant energy 
condition given by Eq. (JTU|) . 

When the transmission coefficient reaches the unit value, that is for E = ?R.E res , 
where E res is the resonance energy, the phase shift S res (E) reaches the value 7r/2. In the 
vicinity of the resonant value, provided that the Breit-Wigner formula [22] is valid, the 
resonance energy is a simple zero of the Jost function f^(E). In the neighborhood of 
E = 3? E res we have 

f+\ E )*[^Pj (E~E res ). (29) 
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Figure 5. Phase shift 6 as a function of the real part of the resonant energy for 
3.145533 < V < 3.15 with a = 1.57 
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If S E res is sufficiently close to the real axis for the approximation (I29p to be applicable 
in a real neighborhood of E = 9ft E res we will have 

8(E) = -arg(U(E)) E=Erea « 6 bg (E) + 8 res (E) (30) 

hence 

8 res (E) « 6(E) - 8 bg (E), (31) 

where 6b g = —org (^H is the background phase shift. Fig [3] shows the 

\ / E=E res 

dependence of the phase shift 6 on the real part of the energy E. For the resonance 
value $IE = —1.0032 we have that 6 = 2.62769, and 6b g = 1.05706, obtaining in this 
way S res = 1.57063 ~ n/2. 

It is worth mentioning that as the lowest-bound energy state reaches the 
supercritical value E = —m, the phase shift 8 takes the value 8 = 7r/2, result that is in 
agreement with the phase shift for the one-dimensional Dirac equation in the presence 
of a symmetric potential well [TJ [8] 

The Wigner time r permits one to estimate how long the bound state lives before 
it degenerates into the negative continuum. From Fig. [4] and the solution to Eq. (|T0|) . 
we can estimate the mean life of the resonant state in terms of the half-width T of the 
Lorentzian curve. For dtE = —1.0032 we have QE = 1.765 x 10 -5 , corresponding to a 
Wigner time of r = -f| ~ ^ = 5.667 x 10 5 expressed in natural units. Resonant states 
make the phase shift surpass the value tt/2 with a Wigner time delay r > Fig. [5] 
shows the behavior of the phase shift 8 as a function of the real part of the resonant 
energy, it can be observed that the phase shift increases surpassing the value 7r/2 and 
after reaching a maximum, it monotonically decreases, showing in this way a behavior 
that has been observed in square well case. [15] 

4. Concluding remarks 

The relation (1101) shows that the one dimensional cusp potential (pQ) supports 
supercritical states and is strong enough to produce resonant states. This result is non 
trivial and interesting in view of the fact that the Dirac equation in a one-dimensional 
vectorial dirac delta interaction does not exhibit resonances [17] . The potential ([1]) does 
not exhibit a square barrier limit, and for very small values of a and constant value 
of 2Va it can be regarded as a delta potential [IB] , therefore the the relation ( TTTT) also 
holds for a delta potential. Eq. <^ represents a local interaction whose support does 
not vanish anywhere and exhibits supercritical resonant states. 

We have computed the phase shift associated with the scattering of a Dirac particle 
by the potential (OQ), and have shown that resonant states make the phase shift surpass 
the value tc/ 2 with a Wigner time-delay r > 0. The asymptotic behavior of the cusp 
potential has permitted us to derive the transmission and reflection amplitudes t(E) 
and r(E) with the help of the Jost solutions. This approach can be applied in the study 
of resonant states in the presence of potentials not allowing exact solutions of the Dirac 
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equation in terms of special functions, this problem will be discussed in a forthcoming 
publication. 
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